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摘 要 : 为 了 找到 一 条 最 短路 径 ， 并 克服 传统 算法 在 路 径 规划 中 不 适合 离散 域 求解 以 及 收 全 速度 慢 等 问题 ， 提 出 一 种 
进 的 狼 群 算法 。 通 过 在 初始 化 阶段 引入 位 置 -次 序 编码 的 方法 ， 研 究 了 离散 域 的 路 径 优 化 ; 同时 在 和 迭代 过 程 中 引入 二 

次 搜索 来 提高 算法 求解 速度 与 精度 ， 以 实现 在 达到 最 大 选 代 次 数 前 出 现 最 优 解 。 结 果 表 明 ， 改 进 的 狼 群 算法 相 比 已 有 

的 算法 求解 精度 更 高 ， 收 化 速度 更 快 ， 更 加 有 效 地 避免 陷入 局 部 最 优 。 可 见 改进 狼 群 算法 可 以 很 好 的 应 用 于 求解 最 优 

路 径 规划 问题 。 
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Improved wolf group algorithm for solving traveling salesman problem 


Huang Haisong, Ren Zhupeng, Wei Jian’an 
(Key laboratory of Advanced Manufacturing Technology of Ministry of Education Guizhou University, Guiyang 550025, China) 


Abstract: In order to find a shortest path and overcome the problems that traditional algorithms are not suitable for discrete 
domain solution and slow convergence in path planning, This Paper proposed an improved wolf group algorithm. By introducing 
the position-order coding method in the initialization phase, this paper studied the path optimization problem of discrete domains. 
At the same time, it introduced a secondary search in the iterative process to improve the speed and accuracy of the algorithm 
to the optimal solution that is before the maximum number of iterations. The results show that the improved wolf group algorithm 
has higher accuracy and faster convergence than the existing algorithms, and it is more effective to avoid falling into local 
optimum. 
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解 的 概率 ， 并 且 缩小 了 解 空 间 。 但 很 多 算法 也 存在 容易 陷入 局 
部 极 值 和 收敛 速度 慢 等 不 足 。 因 此 ， 路 径 规划 问题 叫 依 然 是 具 
有 难度 跟 研 究 价值 的 问题 之 一 ， 

于 狼 群 算法 ， 国 内 外 很 多 学 者 提出 了 不 同 的 改进 方法 。 
献 [8] 通 过 设计 新 的 猛 狼 位 置 更 新 公式 ， 优 化 传统 狼 群 算法 步 
的 方法 进行 改进 。 文献 [9] 通 过 引入 差分 进化 中 的 交叉 、 变 异 、 


旅行 商 问题 (traveling salesman problem，TSP) 是 一 种 常 
见 的 路 径 优 化 问题 ， 其 目的 是 为 了 求 得 一 条 经 过 所 有 城市 的 最 
短路 径 。 现 实生 活 中 , 很 多 问题 都 被 抽象 为 TSP 问题 进行 求解 。 
例如 机 器 人 控制 扬 、 车 间 调 度 、 无 人 机 航 迹 规划 馈 、 计 算 机 联 
网 、 网 络 路 由 器 布设 等 。 虽 然 TSP 属于 NP-hard 册 的 组 合 优化 择 等 操作 ， 提 高 了 算法 的 寻 优 性 能 。 但 是 他 们 都 存在 这 收敛 
问题 ,求解 困难 , 但 是 其 研究 具有 实际 价值 和 重要 的 理论 意义 。 度 慢 ， 易 陷入 局 部 最 优 的 缺点 ， 尤 其 是 未 考虑 到 狼 群 算法 在 
里 论 上 ， 较 小 规模 的 TSP 可 以 用 传统 算法 获得 最 优 解 ， 较 加 和 散 域 求 解 的 问题 。 因 此 为 实现 离散 域 的 路 径 优 化 ,采用 一 种 引 
大 规模 的 TSP 一 般 采 用 启发 式 算法 ,但 启发 算法 也 存在 着 有 限 入 位 置 -次 序 编码 的 方法 。 为 了 提高 狼 群 算法 求解 的 速度 与 精度 ， 
时 间 内 无 法 取得 最 优 解 的 局 限 。 而 由 自然 界 生 物 群 体 演化 而 来 在 迭代 过 程 中 采用 引入 随机 单 点 插入 算 子 、 多 点 插入 算 子 的 方 
的 群体 智能 算法 〈 遗 传 算法 、 粒 子 群 优 化 算法 、 蚁 群 算法 器 、 法 ， 通 过 二 次 搜索 来 实现 在 达到 最 大 迭代 次 数 时 出 现 最 优 解 。 
人 工蜂 群 算法 、 莹 火 虫 群 优化 算法 、 生 物 地 理 迁 移 算 法 、 帝 国 、 
主义 竞争 算法 等 ) ， 为 求解 复杂 路 径 优化 问题 提供 了 不 同 的 方 1 狼 属 算法 的 概述 
法 。 新 兴 的 群体 智能 仿生 算法 名 提高 了 在 有 限时 间 内 获得 最 优 探 狼 游 走 、 召 唤 头 狼 、 狐 狼 围攻 这 3 种 智能 行为 构成 了 狼 
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群 算法 外 的 主体 。 其 规则 是 “优胜 劣 汰 ”的 狼 群 更 新 规则 以 及 
“ 胜 者 为 王 ” 的 头 狼 角逐 规则 。 其 步骤 如 下 ; 
a) 狼 群 的 空间 坐标 在 解 空间 中 随机 初始 化 , 头 狼 即 为 最 优 
标 函 数值 的 人 工 狼 。 
b) 在 搜索 猎物 时 , 探 狼 开始 游 走 。 如 果 探 狼 发 现 的 气味 浓度 
比 头 狼 发 现 的 浓度 大 ， 则 探 狼 处 于 主导 地 位 ， 同 时 更 新 头 狼 位 
。 反 之 ， 探 狼 需要 继续 游 走 ， 直 到 发 现 的 猎物 气味 浓度 大 于 
头 狼 或 者 是 达到 最 大 的 游 走 次 数 ， 那 么 狼 头 将 在 现 有 的 位 置 进 
行 召 唤 。 
9 在 头 狼 发 出 召唤 以 后 , 猛 狼 开始 以 最 大 步 长 向 头 狼 移 动 。 
在 这 个 过 程 中 ， 如 果 发 现 猛 狼 的 目标 值 优 于 头 狼 ， 那 么 将 更 新 
头 狼 。 反 之 ， 猛 狼 将 不 断 移动 ， 直 到 进入 围攻 范围 。 
qd) 联合 探 狼 以 后 ,靠近 头 狼 的 猛 狼 将 围捕 猎物 (把 头 狼 位 置 
视 为 猎物 )。 如 果 在 实施 围攻 以 后 , 发 现 人 工 狼 的 目标 值 优 于 头 
狼 ， 则 更 新 头 狼 位 置 ， 继 续 进 行 捕获 。 反 之 , 则 不 变 。 
e 将 目标 函数 值 较 小 的 人 工 狼 进行 淘汰 ， 并 且 新 的 人 工 狼 
将 在 解 空间 中 随机 生成 ， 不 断 实现 狼 群 的 更 新 。 
人 最 后 头 狼 的 目标 值 需要 根据 是 否 收敛 于 全 局 最 优 解 或 者 
迭代 次 数 是 否 达到 最 大 来 判断 。 如 果 没 有 达到 ， 则 需要 继续 迁 
代 ， 直 到 满足 条 件 。 在 达到 要 求 以 后 ， 输 出 的 目标 函数 值 就 是 
最 优 解 ， 被 求 函数 的 最 佳 位 置 即 为 头 狼 的 空间 坐标 外。 


2 ”改进 狼 群 算法 求解 路 径 规划 


狼 群 算法 虽然 可 以 处 理 简单 的 路 径 规划 问题 ， 但 不 能 实现 
离散 域 的 路 径 规 划 ， 且 和 迭代 速度 较 慢 。 因 此 提出 一 种 改进 狼 群 
算法 , 该 算法 通过 引入 位 置 -次 序 编码 的 方法 ,来 实现 离散 域 的 
路 径 优化 ， 通 过 二 次 搜索 来 提高 狼 群 算法 求解 的 速度 与 精度 ， 
具体 过 程 如 下 。 

2.1 改进 狼 群 算法 求解 路 径 规划 原理 
狼 群 算法 与 路 径 规 划 问 题 的 对 应 关系 如 下 : 
a) 狼 群 算法 一 一 路 径 规划 ; 


J 


Wolve(n).location(k +1) = Wolve(n).location(k) 
+Wolve(n).location(k)*(SDR*(rand*2 —1)(1) 
在 跟踪 模式 下 , 狼 个 体 的 每 次 迭代 都 将 处 于 全 局 最 优 状 态 ， 
如 此 不 断 有 逼近 全 局 最 优 解 .在 此 进程 中 ,每 只 狼 的 速度 标记 为 


人 只, 狼 个 体 分 别 利 


式 (2) (3) 来 更 新 自己 的 速 


度 与 位 置 。 


VO+1)=v ?+e*rand(sd (1)—s0 0 (2) 


其 中 : vO (r+D) ,表示 第 K 只 狠 第 d 位 的 速度 ， 且 d=1,2…,1,，c 


为 常数 ; s(o) (7) 为 狼 群 全 局 最 优 解 ;参考 文献 [10],c=2;rand 为 [0,1] 
的 随机 数 。 
sO(1+1)=s ?+s (+D) G3) 


其 中 ， sD) 表示 当前 位 置 ;更 新 后 的 位 置 是 O441) 。 


记录 当前 迭代 的 最 优 解 并 且 计 算 狼 个 体 的 适应 度 ,继续 进 
行规 定 的 迭代 次 数 。 
2.2 ”改进 狼 群 算法 的 描述 
2.2.1 狼 的 位 置 -次 序 编码 过 程 

在 连续 域 的 函数 优化 中 WPA 可 以 应 用 ,但 是 却 不 适合 在 离 
散 域 TSP 求解 中 直接 应 用 00。 主 要 原因 是 TSP 往往 是 由 自然 
数 表 示 ， 但 搜索 模式 下 变异 发 生 在 实数 域 范围 ， 为 了 克服 此 问 
-次 序 编码 。 


题 ， 引 进位 


狼 群 算法 的 第 t+ 代 第 i 个 体 狼 如 式 (4) 所 示 ，ord( X;(D)) 


表示 将 5,(7),s,,(D),…，s,,(D) 按 升 序 排列 后 w(D 在 所 得 序列 中 的 
次 序 ， 则 式 (5) 是 其 对 应 的 可 行路 径 编 码 ， 这 个 变换 过 程 就 是 
位 置 -次 序 编码 。 


$;(1) = si(D), ;a(t),....,5;, (D (4) 
Wt = (ord(s,,(t)),ord(s,,(t)),...,ord(s,, (t)) (5) 
例如 ，$(b = (2.0337,-19720,4.9901,-4.0015 则 由 上 述 的 


bp) 狼 个 体 一 一 可 行路 径 ; 

9 狼 位 置 变动 一 路径 顺序 变动 ; 

d) 狼 适应 度 一 一 路 径 长 度 。 

在 搜索 模式 时 ， 狼 群 定义 了 三 个 参数 : 个 体 改 变 的 基因 的 
个 数 CDC; 放置 变异 状态 的 记忆 池 SMP; 个 体 上 每 个 基因 芯 


的 
改变 范围 SDR. 具 体 的 搜索 模式 是 ， a) 狼 群 中 当前 状态 的 狼 个 
体 置 于 SMP 并 复制 ， b) 将 变异 发 生 后 的 记忆 池 中 复制 的 副本 
改变 ,并 更 新 到 新 的 位 置 ，c〉 计算 副本 的 适应 度 值 , 狼 最 终 的 更 
新 状态 就 是 其 中 的 一 个 最 优 的 状态 。 

设 狼 个 体 i 的 当前 位 置 为 首先 对 于 狼 群 个 体 
:0O-enowaO ws OO 在 记忆 池 中 复制 其 位 置 然后 在 狼 群 搜 
索 模式 下 选择 SMP-1 个 位 置 进行 变异 操作 , 其 中 每 个 个 体 选择 
CDC 个 基因 位 置 ,变异 公式 如 式 (1) 所 示 。 最 后 对 应 的 可 行路 
径 编码 即 为 SMP 个 位 置 转换 而 成 ， 并 计算 适应 度 ， 更 新 个 体 。 
在 该 搜索 方式 下 所 有 狼 个 体会 逐渐 收敛 于 全 局 最 优 解 。 


WT 


位 置 -次 序 编码 变换 过 程 可 知 ， 对 应 的 位 置 一 次 序 编码 向 量 为 
W(t)=(2,3,L4) 。 


2.2.2 适应 度 的 确定 

设 pi,py,.….，p, (其 中 1,2,...n 为 城市 编号 ) 为 一 条 可 行 访问 
路 径 ， 则 式 (6) 是 路 径 的 总 长 度 。 在 狼 群 算法 中 ,评估 个 体 是 根 
据 路 径 的 长 度 ，d 越 大 ,表示 个 体 适 应 度 越 差 ， 即 路 径 长 度 的 倒 
数 是 个 体 适应 度 (fitness=1/d)。 


7 一 | 
d= APPitD) + dp PD) (6) 
i 


其 中 ，d(p,p,) 一 表示 p, 与 p, 之 间 的 距离 。 
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2.2.3 二 次 搜索 策略 

搜索 模式 是 狼 群 算法 的 主要 优化 过 程 ， 但 它 同 时 也 存在 着 
一 些 问题 ,比如 变异 策略 比较 单一 、 收 敛 性 达 不 到 最 佳 等 方面 。 
当 我 们 要 优化 的 种 群 数目 足够 大 时 ， 因 其 最 大 迭代 次 数 有 限 会 
使 得 优化 变异 后 的 路 径 很 难 出 现 较 优 解 。 当 路 径 没 有 获得 更 优 
解 时 ， 在 搜索 模式 优化 完成 后 ， 为 了 更 好 地 解决 TSP， 引 入 变 
异 算 子 不 同 但 与 狼 群 搜索 模式 相似 的 二 次 搜索 优化 策略 。 如 图 


1 所 示 。 
Y 


最 大 和 迭代 数 ? 


图 1 改进 狼 群 算法 求解 路 径 规 划 的 流程 图 
二 次 搜索 优化 策略 的 两 种 方式 如 下 : 

1) 随机 单 点 插入 算 子 的 描述 

编码 所 有 的 个 体 狼 , 设 c,,c，( 计 为 随机 选择 的 两 个 基因 位 
) 为 其 位 上 的 城市 。 在 


Ad¢-d(c,cn)+d(c,cjw)>d(c,c))+d(c,scn) 式 中 ， 两 座 城 


市 (a,c, ) 的 距离 为 d(ei,c,) 。 为 了 使 整个 路 径 的 距离 得 到 改善 ， 
当 Ad>0， 就 将 j 与 i41 两 基因 位 互 换 。 这 种 方法 任意 选取 其 
中 的 两 个 位 点 进行 基因 交换 ， 对 于 “模式 转化 ”影响 较 小 ， 计 
算 也 并 不 复杂 。 通 过 此 方法 ， 在 每 一 次 迭代 时 都 可 以 选择 距离 
自己 最 近 且 路 径 最 优 的 点 。 如 图 2 所 示 


插入 前 


插入 后 [io[s14[9[1[6[513[7[12?| 


图 2 随机 单 点 插入 
2) 随机 多 点 插入 算 子 的 描述 
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图 3 随机 多 点 插入 
此 种 插入 方法 优点 为 : 
a) 可 以 通过 选择 距离 最 短 的 点 来 改善 算法 的 适应 有 
具有 局 部 优化 作用 ; 
b) 可 以 通过 随机 选取 插入 点 ,来 丰富 可 行 解 的 多 样 性 ; 
0o) 通 过 整体 插入 来 避免 破坏 已 经 存在 的 全 局 最 优 的 部 分 路 
径 序 列 。 


3 ”改进 狼 群 算法 的 仿真 实验 及 分 析 


选取 TSPLIB 标准 库 中 的 几 个 实例 进行 实验 ， 验 证 离散 狼 
群 算 法 的 准确 性 、 有 效 性 和 可 行 性 。 选 择 参考 文献 [14~16] 中 优 
化 算法 的 结果 进行 比较 。 

3.1 实验 结果 
3.1.1 Dantzig42 问题 

针对 Dantzig42 问题 ,本 文 算法 优化 路 径 如 图 4 所 示 ，10 次 

优化 结果 如 表 第 1 列 所 示 。 


bE 
注 


0 20 和 60 80 100 120 140 160 180 


图 4 ” Dantzig42 问题 优化 路 径 


3.1.2 Eil51 问题 
针对 Eil51 问题 ， 本 文 算 法 优化 优化 路 径 如 图 5 所 示 ，10 
次 优化 结果 如 表 1 第 2 列 所 示 。 


$0 


0 10 20 30 和 0 60 70 


图 5 Fil51 问题 优化 路 径 


3.1.3 Berlin52 问题 
针对 Berlin52 问题 ， 本 文 算法 优优 化 路 径 如 图 6 所 示 ， 黑 
色 圆 圈 为 局 部 放大 部 分 ,证 明 没有 出 现 交 叉 。10 次 优化 结果 如 


编码 所 有 的 狼 个 体 ， 不 妨 设 k 为 城市 位 置 的 编号 。 基 因 位 
i 是 随机 选取 的 ， 距 离 k 城 市 最 近 的 t 城 市 是 在 基因 位 以 后 的 n 


一 i 个 位 置 。 其 基因 位 为 j， 将 j 基因 位 以 后 的 n-j 位 整体 插入 


到 基因 位 i 后 ,这 种 线段 算 子 的 插入 , 可 以 大 大 缩减 计算 时 间 ， 
以 最 小 的 迭代 次 数 找到 最 优 解 。 如 图 3 所 示 。 
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表 1 第 3 列 所 示 。 
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图 6 ”Berlin52 问题 优化 路 径 


3.2 ”实验 分 析 
图 4~6 可 得 ,上 述 三 个 TSPLIB 实例 经 过 本 文 IWPA 算法 
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